Open-field PET: Simultaneous brain functional imaging and behavioural response measurements in freely moving small animals.

A comprehensive understanding of how the brain responds to a changing environment requires techniques capable of recording functional outputs at the whole-brain level in response to external stimuli. Positron emission tomography (PET) is an exquisitely sensitive technique for imaging brain function but the need for anaesthesia to avoid motion artefacts precludes concurrent behavioural response studies. Here, we report a technique that combines motion-compensated PET with a robotically-controlled animal enclosure to enable simultaneous brain imaging and behavioural recordings in unrestrained small animals. The technique was used to measure in vivo displacement of [11C]raclopride from dopamine D2 receptors (D2R) concurrently with changes in the behaviour of awake, freely moving rats following administration of unlabelled raclopride or amphetamine. The timing and magnitude of [11C]raclopride displacement from D2R were reliably estimated and, in the case of amphetamine, these changes coincided with a marked increase in stereotyped behaviours and hyper-locomotion. The technique, therefore, allows simultaneous measurement of changes in brain function and behavioural responses to external stimuli in conscious unrestrained animals, giving rise to important applications in behavioural neuroscience.

Direct Estimation of Voxel-Wise Neurotransmitter Response Maps From Dynamic PET Data.

Computational methods, such as the linear parametric neurotransmitter PET (lp-ntPET) method, have been developed to characterize the transient changes in radiotracer kinetics in the target tissue during endogenous neurotransmitter release. In this paper, we describe and evaluate a parametric reconstruction algorithm that uses an expectation maximization framework, along with the lp-ntPET model, to estimate the endogenous neurotransmitter response to stimuli directly from the measured PET data. Computer simulations showed that the proposed direct reconstruction method offers improved accuracy and precision for the estimated timing parameters of the neurotransmitter response at the voxel level ( td=1±2 min, for activation onset bias and standard deviation) compared with conventional post reconstruction modeling ( td=4±7 min). In addition, we applied the proposed direct parameter estimation methodology to a [11C]raclopride displacement study of an awake rat and generated parametric maps illustrating the magnitude of ligand displacement from striatum. Although the estimated parametric maps of activation magnitude obtained from both direct and post reconstruction methodologies suffered from false positive activations, the proposed direct reconstruction framework offered more reliable parametric maps when the activation onset parameter was constrained.

Determining Glucose Metabolism Kinetics Using 18F-FDG Micro-PET/CT.

This paper describes the use of 18F-FDG and micro-PET/CT imaging to determine in vivo glucose metabolism kinetics in mice (and is transferable to rats). Impaired uptake and metabolism of glucose in multiple organ systems due to insulin resistance is a hallmark of type 2 diabetes. The ability of this technique to extract an image-derived input function from the vena cava using an iterative deconvolution method eliminates the requirement of the collection of arterial blood samples. Fitting of tissue and vena cava time activity curves to a two-tissue, three compartment model permits the estimation of kinetic micro-parameters related to the 18F-FDG uptake from the plasma to the intracellular space, the rate of transport from intracellular space to plasma and the rate of 18F-FDG phosphorylation. This methodology allows for multiple measures of glucose uptake and metabolism kinetics in the context of longitudinal studies and also provides insights into the efficacy of therapeutic interventions.

Motion compensation using origin ensembles in awake small animal positron emission tomography.

In emission tomographic imaging, the stochastic origin ensembles algorithm provides unique information regarding the detected counts given the measured data. Precision in both voxel and region-wise parameters may be determined for a single data set based on the posterior distribution of the count density allowing uncertainty estimates to be allocated to quantitative measures. Uncertainty estimates are of particular importance in awake animal neurological and behavioral studies for which head motion, unique for each acquired data set, perturbs the measured data. Motion compensation can be conducted when rigid head pose is measured during the scan. However, errors in pose measurements used for compensation can degrade the data and hence quantitative outcomes. In this investigation motion compensation and detector resolution models were incorporated into the basic origin ensembles algorithm and an efficient approach to computation was developed. The approach was validated against maximum liklihood-expectation maximisation and tested using simulated data. The resultant algorithm was then used to analyse quantitative uncertainty in regional activity estimates arising from changes in pose measurement precision. Finally, the posterior covariance acquired from a single data set was used to describe correlations between regions of interest providing information about pose measurement precision that may be useful in system analysis and design. The investigation demonstrates the use of origin ensembles as a powerful framework for evaluating statistical uncertainty of voxel and regional estimates. While in this investigation rigid motion was considered in the context of awake animal PET, the extension to arbitrary motion may provide clinical utility where respiratory or cardiac motion perturb the measured data.

List-mode image reconstruction for positron emission tomography using tetrahedral voxels.

Image space decomposition based on tetrahedral voxels are interesting candidates for use in emission tomography. Tetrahedral voxels provide many of the advantages of point clouds with irregular spacing, such as being intrinsically multi-resolution, yet they also serve as a volumetric partition of the image space and so are comparable to more standard cubic voxels. Additionally, non-rigid displacement fields can be applied to the tetrahedral mesh in a straight-forward manner. So far studies incorporating tetrahedral decomposition of the image space have concentrated on pre-calculated, node-based, system matrix elements which reduces the flexibility of the tetrahedral approach and the capacity to accurately define regions of interest. Here, a list-mode on-the-fly calculation of the system matrix elements is described using a tetrahedral decomposition of the image space and volumetric elements-voxels. The algorithm is demonstrated in the context of awake animal PET which may require both rigid and non-rigid motion compensation, as well as quantification within small regions of the brain. This approach allows accurate, event based, motion compensation including non-rigid deformations.

Full field spatially-variant image-based resolution modelling reconstruction for the HRRT.

Accurate characterisation of the scanner's point spread function across the entire field of view (FOV) is crucial in order to account for spatially dependent factors that degrade the resolution of the reconstructed images. The HRRT users' community resolution modelling reconstruction software includes a shift-invariant resolution kernel, which leads to transaxially non-uniform resolution in the reconstructed images. Unlike previous work to date in this field, this work is the first to model the spatially variant resolution across the entire FOV of the HRRT, which is the highest resolution human brain PET scanner in the world. In this paper we developed a spatially variant image-based resolution modelling reconstruction dedicated to the HRRT, using an experimentally measured shift-variant resolution kernel. Previously, the system response was measured and characterised in detail across the entire FOV of the HRRT, using a printed point source array. The newly developed resolution modelling reconstruction was applied on measured phantom, as well as clinical data and was compared against the HRRT users' community resolution modelling reconstruction, which is currently in use. Results demonstrated improvements both in contrast and resolution recovery, particularly for regions close to the edges of the FOV, with almost uniform resolution recovery across the entire transverse FOV. In addition, because the newly measured resolution kernel is slightly broader with wider tails, compared to the deliberately conservative kernel employed in the HRRT users' community software, the reconstructed images appear to have not only improved contrast recovery (up to 20% for small regions), but also better noise characteristics.

Impact of extraneous mispositioned events on motion-corrected brain SPECT images of freely moving animals.

PURPOSE: Single photon emission computed tomography (SPECT) brain imaging of freely moving small animals would allow a wide range of important neurological processes and behaviors to be studied, which are normally inhibited by anesthetic drugs or precluded due to the animal being restrained. While rigid body motion of the head can be tracked and accounted for in the reconstruction, activity in the torso may confound brain measurements, especially since motion of the torso is more complex (i.e., nonrigid) and not well correlated with that of the head. The authors investigated the impact of mispositioned events and attenuation due to the torso on the accuracy of motion corrected brain images of freely moving mice. METHODS: Monte Carlo simulations of a realistic voxelized mouse phantom and a dual compartment phantom were performed. Each phantom comprised a target and an extraneous compartment which were able to move independently of each other. Motion correction was performed based on the known motion of the target compartment only. Two SPECT camera geometries were investigated: a rotating single head detector and a stationary full ring detector. The effects of motion, detector geometry, and energy of the emitted photons (hence, attenuation) on bias and noise in reconstructed brain regions were evaluated. RESULTS: The authors observed two main sources of bias: (a) motion-related inconsistencies in the projection data and (b) the mismatch between attenuation and emission. Both effects are caused by the assumption that the orientation of the torso is difficult to track and model, and therefore cannot be conveniently corrected for. The motion induced bias in some regions was up to 12% when no attenuation effects were considered, while it reached 40% when also combined with attenuation related inconsistencies. The detector geometry (i.e., rotating vs full ring) has a big impact on the accuracy of the reconstructed images, with the full ring detector being more advantageous. CONCLUSIONS: Motion-induced inconsistencies in the projection data and attenuation/emission mismatch are the two main causes of bias in reconstructed brain images when there is complex motion. It appears that these two factors have a synergistic effect on the qualitative and quantitative accuracy of the reconstructed images.

Evaluation of a direct 4D reconstruction method using generalised linear least squares for estimating nonlinear micro-parametric maps.

OBJECTIVE: Estimation of nonlinear micro-parameters is a computationally demanding and fairly challenging process, since it involves the use of rather slow iterative nonlinear fitting algorithms and it often results in very noisy voxel-wise parametric maps. Direct reconstruction algorithms can provide parametric maps with reduced variance, but usually the overall reconstruction is impractically time consuming with common nonlinear fitting algorithms. METHODS: In this work we employed a recently proposed direct parametric image reconstruction algorithm to estimate the parametric maps of all micro-parameters of a two-tissue compartment model, used to describe the kinetics of [[Formula: see text]F]FDG. The algorithm decouples the tomographic and the kinetic modelling problems, allowing the use of previously developed post-reconstruction methods, such as the generalised linear least squares (GLLS) algorithm. RESULTS: Results on both clinical and simulated data showed that the proposed direct reconstruction method provides considerable quantitative and qualitative improvements for all micro-parameters compared to the conventional post-reconstruction fitting method. Additionally, region-wise comparison of all parametric maps against the well-established filtered back projection followed by post-reconstruction non-linear fitting, as well as the direct Patlak method, showed substantial quantitative agreement in all regions. CONCLUSIONS: The proposed direct parametric reconstruction algorithm is a promising approach towards the estimation of all individual microparameters of any compartment model. In addition, due to the linearised nature of the GLLS algorithm, the fitting step can be very efficiently implemented and, therefore, it does not considerably affect the overall reconstruction time.

Isotope specific resolution recovery image reconstruction in high resolution PET imaging.

PURPOSE: Measuring and incorporating a scanner-specific point spread function (PSF) within image reconstruction has been shown to improve spatial resolution in PET. However, due to the short half-life of clinically used isotopes, other long-lived isotopes not used in clinical practice are used to perform the PSF measurements. As such, non-optimal PSF models that do not correspond to those needed for the data to be reconstructed are used within resolution modeling (RM) image reconstruction, usually underestimating the true PSF owing to the difference in positron range. In high resolution brain and preclinical imaging, this effect is of particular importance since the PSFs become more positron range limited and isotope-specific PSFs can help maximize the performance benefit from using resolution recovery image reconstruction algorithms. METHODS: In this work, the authors used a printing technique to simultaneously measure multiple point sources on the High Resolution Research Tomograph (HRRT), and the authors demonstrated the feasibility of deriving isotope-dependent system matrices from fluorine-18 and carbon-11 point sources. Furthermore, the authors evaluated the impact of incorporating them within RM image reconstruction, using carbon-11 phantom and clinical datasets on the HRRT. RESULTS: The results obtained using these two isotopes illustrate that even small differences in positron range can result in different PSF maps, leading to further improvements in contrast recovery when used in image reconstruction. The difference is more pronounced in the centre of the field-of-view where the full width at half maximum (FWHM) from the positron range has a larger contribution to the overall FWHM compared to the edge where the parallax error dominates the overall FWHM. CONCLUSIONS: Based on the proposed methodology, measured isotope-specific and spatially variant PSFs can be reliably derived and used for improved spatial resolution and variance performance in resolution recovery image reconstruction. The benefits are expected to be more substantial for more energetic positron emitting isotopes such as Oxygen-15 and Rubidium-82.

Convergence optimization of parametric MLEM reconstruction for estimation of Patlak plot parameters.

In dynamic positron emission tomography data many researchers have attempted to exploit kinetic models within reconstruction such that parametric images are estimated directly from measurements. This work studies a direct parametric maximum likelihood expectation maximization algorithm applied to [(18)F]DOPA data using reference-tissue input function. We use a modified version for direct reconstruction with a gradually descending scheme of subsets (i.e. 18-6-1) initialized with the FBP parametric image for faster convergence and higher accuracy. The results compared with analytic reconstructions show quantitative robustness (i.e. minimal bias) and clinical reproducibility within six human acquisitions in the region of clinical interest. Bland-Altman plots for all the studies showed sufficient quantitative agreement between the direct reconstructed parametric maps and the indirect FBP (--0.035x+0.48E--5).